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A simple tight-binding model is constructed for the description of the electronic structure of 
some Ce-based filled skutterudite compounds showing an energy gap or pseudogap behavior. 
Assuming band-diagonal electron interactions on this tight-binding model, the optical conduc- 
tivity spectrum is calculated by applying the second-order self-consistent perturbation theory 
to treat the electron correlation. The correlation effect is found to be of great importance on 
the description of the temperature dependence of the optical conductivity. The rapid disap- 
pearance of an optical gap with increasing temperature is obtained as observed in the optical 
experiment for Ce-based filled-skutterudite compounds. 
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■ 1. Introduction 



In recent years, the filled skutterudite compounds with 
chemical formula RT4X12 (R = rare-earth elements, T = 
Fe, Ru and Os, and X = P, As, and Sb) have attracted 
much attention because of several interesting physical 
properties including a new type of heavy-fermion super- 
conductivity observed in Pr-based compounds. 1 In these 
filled skutterudites, some Ce-based compounds have re- 
ceived interest in relation to the Kondo insulators. 2 
For CeRu4Sbi2, for example, it was reported that the 
(pseudo-)gap (~ lOmeV) opens in the spectra obtained 
by several measurements: photoemission, 3 optical con- 
ductivity 4 and inelastic neutron scattering, 5,6 though the 
temperature dependence of resistivity shows metallic be- 
havior. 7 For CeOs4Sbi2, a semiconducting behavior in 
the temperature dependence of resistivity was reported, 8 
and a gap formation (E g ~ lOmeV) was observed in 
the optical measurement 9 at low temperatures. In ad- 
dition, a large peak at 70 meV, which is called a mid- 
infrared (MIR) peak, was also observed. With increas- 
ing temperature, the gap is gradually filled up. At high 
temperatures, the peak structure loses its intensity and 
the gap disappears. 9 The overall structure of the opti- 
cal conductivity and the temperature dependence of the 
spectrum are very similar to those of CeRu4Sbi2 (MIR 
peak at O.leV) in spite of the difference in the transport 
properties. 4 It should be noted that these overall struc- 
tures are also found in a typical Kondo insulator YbBi2 10 
{E g ~ 20meV and MIR peak at 0.3 eV). 

In fact, the overall band structure obtained by the 
band calculation for CeOs4Sbi 2 11 is qualitatively simi- 
lar to that for CeRu4Sbi 2 , 12 though there is a discrep- 
ancy that the gap opens in CeRu4Sbi 2 which shows 
metallic temperature dependence in transport proper- 
ties, whereas it does not open in CeOs4Sbi 2 whose re- 
sistivity shows a semiconducting behavior. This discrep- 
ancy is not yet resolved and the difference in the trans- 
port properties might originate from unknown mecha- 
nism at low energies. They might be attributed to im- 
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purities or local strain. In that case, refinement of the 
sample quality may resolve the discrepancy. Neverthe- 
less, one can expect that similar structures of optical 
conductivity spectra for both Ce-based filled skutteru- 
dite compounds originate from similar band structures. 

In this paper, we construct a simple tight-binding 
model which reproduces the overall structure of the 
band dispersion obtained by the band calculation for 
CeRu4Sbi2, 12 and introduce the electron-electron inter- 
actions in the tight-binding model. We investigate the 
electron correlation effect on the temperature depen- 
dence of optical conductivity spectra and try to compare 
with the experiments. 

2. Model and Formulation 

Excluding T sites in RT4X12, the crystal structure of 
the filled skutterudites is regarded as the body-centered 
cubic (bec) structure consisting of X12 icosahedrons 
with a rare-earth ion R at a center. In CeOs4Sbi2 or 
CeRu4Sbi2, the f electron orbitals of Ce are expected 
to hybridize with the f symmetry combination of the 
p orbitals of each Sbi2 cluster. Assuming an effective 
overlap integral between f wave functions on nearby R- 
sites through Sbi2 clusters, we consider a tight-binding 
model (with nearest-neighbor (ff(r)-bonds) 13 consisting 
of f symmetry wave functions on bec lattice sites. By the 
spin-orbit interaction, the 14-fold f states are split into 
the total angular momentum J = 5/2 and J — 7/2 eigen- 
states. Besides, under Th symmetry, the J = 5/2 states 
are split into T 5 and r 6 7 states, and the J = 7/2 states 
are split into two Ts's and T67. 14 Although overlap inte- 
grals (Slater-Koster parameters) for f electrons have been 
obtained in ref. 15 without the spin-orbit interaction, we 
need an expression in terms of T a states for each J value. 
This scheme was performed in ref. 16 for YbBi2 and we 
follow the same procedure here. The resulting f bands, 
which will be shown in the next section, are rather sim- 
ilar to those above the Fermi level obtained by the LDA 
band calculation. 12 

The LDA bands below the Fermi level are rather com- 
plicated. We found, however, that the top-most band (to- 
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Fig. 1. Schematic picture of the energy gap structure, (a) The 
present model, (b) The typical situation of the Kondo insulator. 



tal width of 0.07eV) has a flat part of the width of about 
0.02eV along T-N-P-r lines in the Brillouin zone, which 
is comparable to that of the f-bands above the Fermi 
level. This may be partly due to the mixing with the 
f states above the gap since it contains a considerable f 
components. Therefore, we express it by the simple tight- 
binding band of f symmetry states composed of Sbi2 p 
orbitals. Note that the spin-orbit interaction can be ne- 
glected in this Sbi2 p band since it is very weak in the 
p band. Under Th group, the f symmetry states with- 
out the spin-orbit interaction are split into A u and two 
T u 's. In this paper, we assume that the Sbi2 p band with 
f A u symmetry lies below the Fermi level and the low- 
energy optical excitation originates from the transition 
from this to f bands above the Fermi level for CeOs4Sbi2 
or CeRu4Sbi2- In analogy with semiconductors, the for- 
mer and the latter correspond to the valence and the con- 
duction band, respectively. It should be noted that the 
gap which we consider in the present model is not the hy- 
bridization gap in contrast to the gap observed in other 
heavy fermion compounds classified into the Kondo in- 
sulator / semiconductor, e.g. YbBi2. 10 The difference of 
the (indirect) energy gap structure between the present 
model and the typical Kondo insulator is schematically 
shown in Fig. 1. In the present model, the 4f band (which 
is schematically shown as an energy level in Fig. 1(a)) sits 
above the valence band, and the gap opens between the 
4f band and the valence band. On the other hand, in the 
typical Kondo insulator, the gap is constructed by the hy- 
bridization between the 4f level and the conduction band 
(Fig. 1(b)). It should be noted that there should be some 
hybridization between the 4f and the valence bands also 
in the case (a). In the present calculation, therefore, we 
assume the weak hybridization between the Sbi2 p band 
and f bands, hence our model lies between the cases (a) 
and (b). 

We consider the effective tight-binding model men- 
tioned above as a free system without electron-electron 
interactions and assume effective Coulomb and exchange 
interactions as follows: 19 



^ — Wband+^int, 

Wband = ^ E ^k C 7kcr C 7kg-i 
7 k,(T 

n int = u^Y, 

7 i 



(i) 

(2) 



7<7' i.a 



t/3 ^ ^ Tl^i a Tl^> id 
7<7' ?,cr 



where E£ denotes the diagonalized tight-binding band 
energy for the band 7. We denote the annihilation (cre- 
ation) operator for the band 7 in the site representation 
as c 7 i CT (cjj, ia .). Here the index a denotes a pair of the 
time-reversal states which we call spin hereafter. We as- 
sume the the para state for the spin and omit the index 
a. In the Hamiltonian, we have introduced the above 
form of electron interactions expressed by the param- 
eters U, U2, U3 and J. They correspond to the intra- 
band Coulomb interaction, the inter-band anti-parallel 
spin Coulomb interaction, the inter-band parallel spin 
Coulomb interaction, and the inter-band exchange inter- 
action, respectively. In the interaction parts Hint of the 
Hamiltonian, we assume band-diagonal electron interac- 
tions. Usually the interactions are introduced between 
the localized electron orbitals before the diagonalization 
of , but other choices of basis functions are also possi- 
ble in general. They are equivalent with each other, but a 
difference may arise when one neglects some of the matrix 
elements. The present form is suitable for the calculation 
of the many-body effect, as will be seen below. Besides, 
in the present paper, the interaction parameters are cho- 
sen to fit the experiments, so that a different choice of 
the interaction form might not change our result. 

In the Hamiltonian (3), the Coulomb interactions are 
introduced in all bands since there is the hybridization 
between the Sbi2 p band and f bands, although the 
strong electron-electron interactions should exist only be- 
tween the f electrons in real systems. However, the equal 
Coulomb interaction may act more strongly on the re- 
gion of high density of states, so that the f-components 
will be more affected than the others. For simplicity and 
in order to reduce the number of parameters, the fol- 
lowing relations are further assumed; U2 = U — J and 
U 3 = U - 2 J. 17 In this case, Hint can be expressed as 

Hint — C 7'i<T' C 7i(T C 7«o' C 7'*o'' 
77'crcr' i 

77'cto"' i.a 

which has a rotational symmetry in both the orbital and 
the spin states. 

In the present model, it is noted that the electron cor- 
relation effect beyond the Hartree-Fock approximation 
does not appear in the ground state at T = 0. In order 
to take account of the correlation effect, we apply the self- 
consistent second-order perturbation theory (SCSOPT) 
to the present model with the local approximation for the 
self-energy part of the Green's function. It is expected 
that the SCSOPT is sufficient to investigate the correla- 
tion effect in the low-energy and low temperature region. 
Using the density of states (DOS) -D 7 (^) for each band 
7, the local Green's function is expressed as follows; 



G 7 (e + W) 



du 



(5) 
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Fig. 2. Second-order perturbation terms of the self energy (e + 
i<5) for the 7 component. The solid line denotes the full local 
Green's function. Dashed lines indicate electron-electron inter- 
actions: U, U2, U3 or J. 



(6) 



where iV is the number of sites and (5 — ^ + . The self 
energy £ 7 (e + 18) consists of the constant Hartree term 
and the second-order perturbation term £ 7 (e + id) . The 
former is neglected here since it can be regarded as al- 
ready included in the band calculation. The latter is cal- 
culated from the full local Green's function as shown in 
Fig. 2. In these diagrams, the solid line denotes the full 
local Green's function G 7 (e + iS) which is determined 
self-consistently together with the self-energy. 

The number density of electrons is denoted by n 7 : 



p 7 (e) 



1 



dep 7 (e)/(e), 



ImG 7 (e + i<S), 



(7) 



(8) 



where /(e) is the Fermi distribution function : /(e) = 
l/(e^( e -* t ) + l) and (3 = iksT)' 1 . The chemical potential 
fj, should be determined so that n to tai(= ^ 7 n 1 ) is equal 
to the given total density of electrons. In the practical 
calculation, we set the chemical potential at T = in 
the middle point between the top of the valence band 
and the bottom of the conduction band, and we use the 
fixed value of the chemical potential at any temperature 
for simplicity. 

According to the linear response theory, the optical 
conductivity is obtained from the current-current corre- 
lation function (eq.(A-l) in Appendix). In deriving the 
formula, we neglect k and 7 dependence of the velocity 
matrix element of the current operator. The transition 
of electrons in the optical absorption is the direct one 
conserving the momentum. We, however, found that the 
formula eq.(A-l) can not reproduce the experiments at 
all because a sharp threshold behavior appears at the 
direct gap edge (see Fig. A-l in Appendix). We there- 
fore assume that the momentum conservation is violated 
in the real system because of imperfections and phonon- 
assisted transitions. By these simplifications, the current- 
current correlation function is reduced to the joint-DOS- 
type form. 18 In the present paper, we use the following 
expression for the optical conductivity by omitting con- 
stant factors: 



cr(w) 



7,7' 



de/o 7 (e)p 7 (e + oj) 



f(e) -/(e + w) 



0) 



There is a problem that the joint-DOS-type form does 
not give the correct spectrum for ui — > at T > 0. In 



the present study, however, we use the joint-DOS-type 
form of the optical conductivity since we treat the system 
with the gap and investigate the overall structure of the 
spectrum. 

3. Results 

Figure 3 shows the band dispersion and the total DOS 
we use in the case without interactions. The bands above 
the Fermi level (indicated by the dashed line) are the f 
bands constructed by the tight-binding model with the 
two-center integral (fie) = 0.005 Ry. In order to repro- 
duce the band dispersions obtained by the LDA band 
calculation, we choose Ep 5 = 0.870 Ry and Er e7 = 0.867 
Ry for J = 5/2. The band structure just above the Fermi 
level corresponds to J — 5/2 and that in the high-energy 
region corresponds to J = 7/2. J = 7/2 states have ma- 
trix elements for the crystalline electric field (CEF) term 
(Og — Oq) in the Tj, group, but we neglected them for 
simplicity since the J = 7/2 states lie at high energies. 
E a), E ( 2) and E Tb7 are chosen as 0.894 Ry, 0.902 Ry 

and 0.889 Ry, respectively. These values are chosen in- 
dependently from the CEF parameters since there might 
be contributions due to hydridization from implicit sur- 
rounding orbitals (not included in the present model ex- 
plicitly). For the valence band corresponding to Sbi2 p 
band with the A u symmetry, the value of (ffcr) is set 
equal to 0.002 Ry and the energy level to 0.857 Ry. The 
hybridization between the Sbi2 p band and f bands is in- 
troduced by (ffcr)', which is set equal to 0.003 Ry. There 
is an indirect gap between T and H points in accord with 
the LDA band. The position of the direct gap differs from 
the r point in the LDA band, but as will be seen later, 
we consider the indirect gap as most important. 

Concerning the valence band, there is a rather flat part 
between N and P points. This is an effect of the hy- 
bridization (ffcr)' between the valence band and the 4f 
band above Ep and yield a peak in the DOS. The total 
valence band width assumed in the present calculation is 
rather narrow than that estimated from the band calcu- 
lation. 11 As mentioned in the previous section, we have 
regarded the present band with the narrow width and the 
above-mentioned peak at 0.855 Ry as corresponding to 
the sharp peak structure at the top of the valence band 
obtained from the band calculation, and the contribu- 
tions from the bands below it are neglected. It may be 
allowed since the main contribution to the low-energy ex- 
citation spectrum of the optical conductivity originates 
from the bands near the gap. 

In the present calculation, we set U = 0.05eV and 
J = 0.03eV. We choose these values to reproduce the 
temperature dependence of the optical conductivity in 
the experiment. Though the DOS does not have any tem- 
perature dependence within the Hartree approximation, 
the strong temperature dependence shows up if the cor- 
relation effect is included. Figure 4(a) shows the temper- 
ature dependence of the total DOS p(e) (it is noted that 
the chemical potential is set to the origin in the figures of 
the total DOS). The sharp peaks lose their intensity and 
their structures become smooth as the temperature in- 
creases. Note that the temperature-dependent spectrum 
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Fig. 3. The total DOS (left) and the band structure (right) for 
the non-interacting system. The Fermi level is indicated by the 
dashed line. 



(a) 




£»r=o.ooo 


(7=0.05 eV 




(eV) 0.010 


7=0.03 eV 


kj 


0.030 



-0.2 0.2 0.4 0.6 0.8 

E (eV) 



(b) 






" /t B r=o.ooo 






. (eV) 0.008 






0.010 






0.012 






0.015 






0.030 











-0.04 -0.02 0.02 0.04 

£ (eV) 



Fig. 4. (a) Temperature dependence of the total DOS p(e). (b) 
Expanded figure of the total DOS near the gap. 

of the DOS cannot be obtained without taking account of 
the correlation effect. The temperature dependence near 
the gap is shown in Fig. 4(b). One can see that the gap 
is filled up with increasing temperature and it turns into 
the pseudogap at the high temperatures comparable to 
the gap size. The peak structures near the gap lose their 
intensity and they shift to the center of the gap. In the 
present model we use, the gap at T = is equal to about 
0.02 eV. Thus it is recognized that the gap is filled in at 
much lower temperature than that corresponding to the 
gap size itself. At the high temperatures corresponding 
to the gap size, the gap structure disappears. 

The temperature dependence of the optical conduc- 
tivity is shown in Fig. 5. The spectrum has a large peak 
near u> ~ 0.1 eV and the (indirect) gap which reflects the 
DOS gap at 0.02eV. With increasing temperature, the in- 
traband contribution (the Drude part) appears and the 
spectrum has finite intensity in the gap region. The large 
peak loses its intensity with increasing temperature. At 
the highest temperature we calculated (k-gT = 0.030 eV, 
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Fig. 5. Temperature dependence of the optical conductivity cal- 
culated by using the SCSOPT. Inset shows the experimental data 
of CeRu4Sbi2. 4 Peaks due to phonons observed in the gap are 
omitted. Note that the logarithmic scale is used for u) axis. 
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Fig. 6. Comparison of the interband contributions of optical con- 
ductivity spectra calculated by using (a) the SCSOPT and (b) 
the Hartree-Fock (HF) approximation. 

which is higher than the room temperature), the gap 
and the peak structures disappear completely. In the op- 
tical experiment for CeRu4Sbi2, 4 the gap is filled up with 
increasing temperature and it disappears completely at 
300K as shown in the inset of Fig. 5. In the experimental 
data, one can see also the decrease of the intensity of the 
MIR peak with increasing temperature. 

In order to see the correlation effect in the tempera- 
ture dependence of the optical conductivity more clearly, 
we show a comparison of spectra calculated by using the 
SCSOPT and the Hartree-Fock (HF) approximation in 
Fig. 6. Note that the figure shows the interband contri- 
bution ointer^) which consists of the part 7 ^ 7' of the 
summation in eq.(9). Obviously the rapid decrease of the 
intensity corresponding to the MIR peak can not be ob- 
tained within the HF approximation. Since the temper- 
ature dependence of the optical conductivity calculated 
by using the HF approximation originates from only the 
Fermi distribution function, the large peak of the spec- 
trum does not show any remarkable change with increas- 
ing temperature. 

4. Summary and Discussion 

We have calculated the density of states and the opti- 
cal conductivity for the simple tight-binding model, in- 
cluding local Coulomb and exchange interactions by SC- 
SOPT; thereby the correlation effect on the temperature 
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dependence of optical conductivity spectra has been in- 
vestigated. The intensity in the gap region grows rapidly 
and the large peak above the gap shifts to the lower 
energy with increasing temperature and the gap dis- 
appears at high temperatures; the gap is filled at the 
lower temperature than that corresponding to the gap 
size itself. In the optical measurement for CeOs4Sbi2 9 
and CeRu4Sbi2, 4 the gap structure disappears at a lower 
temperature than the gap size. The rapid disappearing of 
the gap and the shift of the peak structure cannot be ob- 
tained without the correlation effect. Thus it is concluded 
that the correlation effect is important to explain the 
temperature dependence of the optical conductivity spec- 
tra obtained by the experiments for the Ce-based filled 
skutterudites : CeOs4Sbi2 and CeRu4Sbi2. In the pho- 
toemission experiment for CeRu4Sbi 2 , 3 the temperature 
dependence of the photoemission spectrum was investi- 
gated. In their result, the pseudo-gap is filled with in- 
creasing temperature. The temperature-dependent DOS 
spectrum cannot be explained without the correlation 
effect as mentioned repeatedly. 

In the present study, it has been shown that the large 
peak corresponding to the MIR peak originates from the 
optical transition to the J = 5/2 states from the va- 
lence band. Thus, there should be the another peak cor- 
responding to the transition to the J = 7/2 states at 
higher energy regions. In fact, in the present result, we 
can obviously see a peak structure corresponding to the 
transition to the J = 7/2 state at about 0.5 eV in Fig. 5. 
When we consider that the large peak (MIR peak) near 
0.1 cV observed in experimental data for CeOs4Sbi 2 and 
CeRu4Sbi2 originates in the optical transition from the 
valence band of Xi 2 clusters to the conduction band con- 
sisting of dispersive f J = 5/2 bands, it can be expected 
that there is another peak structure corresponding to 
the transition to f J = 7/2 bands. In the experimental 
result for CeRu4Sbi2, 4 there is a peak around 1 eV (in- 
set of Fig. 5). For CeOs4Sbi 2 , 9 there seems to be also a 
peak structure near 1 eV. According to the band calcu- 
lation, 11 the magnitude of the spin-orbit splitting can be 
estimated as 0.02 ~ 0.03 Ry, i.e. a few hundreds meV. 
Thus the peak position near leV in the experimental 
spectrum is too high and one can not simply conclude 
that it originates from the transition to J = 7/2 band. 
Although this problem cannot be readily explained by 
the present study, it should be reconciled whether there 
is a possibility that the transition to J = 7/2 corresponds 
to the peak structure near 1 cV in the optical conduc- 
tivity spectrum for CcOs4Sbi2 or CeRu4Sbi2 both from 
experimental and theoretical sides. In order to discuss 
the detail structure of the spectrum, however, we need 
a further study with a more elaborated band-structure 
model. 

Finally, we comment that a similar theoretical scheme 
adopted in the present paper was recently applied to the 
optical spectra of the typical Kondo insulator YbBi 2 20 
with a lot of common features, indicating an underlying 
shared structures in both materials. 
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Appendix: Comment on the joint-DOS-type op- 
tical conductivity 

In this Appendix, we comment on the comparison of 
the optical conductivity spectrum obtained from the lo- 
cal DOS (eq.(9)) and that obtained from the direct opti- 
cal transition. When the momentum conservation is as- 
sumed in derivation of the optical conductivity, we ob- 
tain the following expression for the optical conductivity 
with the direct optical transition (assuming velocity ma- 
trix elements to be constant and neglecting the vertex 
correction again): 



^direct M 



- to/ 



d£ Pk( £ Vk (e + ">) 



f(s)-f(e + u) 



id 



AW 



7r Im e + W-^-ET(e + W)' 



(AT) 



(A-2) 



where we omit the constant factor in defining cr ( jj rect 
In contrast to ddiroct(^), we denote here the joint-DOS- 
type optical conductivity expressed by eq.(9) as <7i OC ai(w). 
The comparison of aiocai(w) and <7di rC ct(w) for the non- 
interacting case is shown in Fig. AT. The gap appearing 
in (7i oca i(u>) reflects the indirect gap and it is less than the 
direct gap reflected in tJdiroct(^), but the overall structure 
including positions of large peaks is consistent with each 
other. 

For T > 0, the gap of DOS is smeared out by the 
many-body effect. In this case, the Drude part for 7 = 
i becomes as crS rudc (w -> 0) = ^i/9 7 (0) from (AT) if 
Pk( £ ) = -i Im e -_E i >+i/(2r) ' althou g h it becomes (p 7 (0)) 2 
from (9). Therefore, it should be noted that the joint- 
DOS-type form of the optical conductivity does not give 
the correct spectrum for ui — > at T > and p 7 (0) > 0. 
Despite this deficiency, oi oca i(w) can explain the gross 
feature of the observed spectra at finite frequencies in 
the realistic situations. 
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